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We present an effective Z 2 gauge theory that captures various competing phases in spin-1/2 
kagome lattice antiferromagnets: the topological Z 2 spin liquid (SL) phase, and the 12-site and 36- 
site valence bond solid (VBS) phases. Our effective theory is a generalization of the recent Z 2 gauge 
theory proposed for SL phases by Wan and Tchernyshyov. In particular, we investigate possible 
VBS phases that arise from vison condensations in the SL. In addition to the 12-site and 36-site 
VBS phases, there exists 6-site VBS that is closely related to the symmetry-breaking valence bond 
modulation patterns observed in the recent density matrix renormalization group simulations. We 
find that our results have remarkable consistency with a previous study using a different Z 2 gauge 
theory. Motivated by the lattice geometry in the recently reported vanadium oxyfluoride kagome 
antiferromagnet, our gauge theory is extended to incorporate lowered symmetry by inequivalent 
up- and down-triangles. We investigate effects of this anisotropy on the 12-site, 36-site, and 6-site 
VBS phases. The 12-site VBS is stable to anisotropy while the 36-site VBS undergoes severe dimer 
melting. Interestingly, any analogue of the 6-site VBS is not found in this approach. We discuss the 
implications of these findings and also compare the results with a different type of Z 2 gauge theory 
used in previous studies. 


I. INTRODUCTION 

Exploring exotic quantum states of matter has been 
one of the main themes in condensed matter physics. 
In the last two decades, frustrated quantum antiferro¬ 
magnets have received intense interests due to the ca¬ 
pability of harboring various exotic spin states. Z 2 spin 
liquids (SL) and valence bond solids (VBS) are promi¬ 
nent examples.The simplest description of the Z 2 
spin liquid can be obtained from a resonating valence 
bond (RVB) wave function which has characteristics of 
Z 2 topological order and fractionalized quasiparticle exci¬ 
tations such as spinons and visons.^^® On the other hand, 
VBS is a crystalline order of valence bonds, and its ele¬ 
mentary excitations are triplons which are confined two- 
spinon states. 

The spin-1/2 nearest-neighbor antiferromagnetic 
Heisenberg model on the kagome lattice has been 
extensively studied as a promising model for exotic spin 
states. The nature of the ground state in this model 
system was under debate due to the existence of various 
competing low energy states. Earlier studies of analytic 
approaches and numerical computations proposed three 
major ground state candidates: the gapped Z 2 
U(l) Dirac and a VBS with a 36-site unit 

cell^®“^^. Recent density matrix renormalization group 
(DMRG) simulations suggest that this model has a 
Z 2 SL ground state characterized by a Z 2 topological 
order and a finite energy gap to excitations.^®’^® The 
simulations also reveal diamond-shaped patterns in 
valence bond correlations of the SL state, indicating that 
the long range ordered counterpart of the state may be a 
VBS with a 12-site unit cell.^^ Although the topological 
Z 2 SL is the ground state of the nearest-neighbor 


Heisenberg model, other states mentioned above may 
also be realized by adding small perturbations to the 
Heisenberg model. For example, the 36-site VBS is 
stabilized by second nearest-neighbor ferromagnetic 
interactions as shown in a previous variational Monte 
Carlo approach.^® With many low energy states in 
competition, the kagome antiferromagnets provide a 
playground for various competing SL and VBS phases. 

On the experimental side, several materials re¬ 
alize the spin-1/2 kagome lattice antiferromag¬ 
net. Herbertsmithite [ZnCu 3 (OH) 6 Cl 2 ],®^^^^ the 
deformed kagome antiferromagnet [Rb 2 Cu 3 SnFi 2 ],^®^®^ 
and the recently discovered vanadium oxyfluoride 
[(NH4)2(C7 Hi4N)(V706Fi 8)]®®’®'* are such examples. In 
these materials, there are generically small perturbations 
to the nearest-neighbor model on the ideal lattice, such 
as further neighbor interactions, Dzyaloshinskii-Moriya 
interactions, and lattice distortions. As an example, 
vanadium oxyfluoride has a kagome lattice formed by 
V^"*" spin-1/2 moments with inequivalent up- and down- 
triangles which lowers the sixfold rotation symmetry 
down to threefold. Such small perturbations can easily 
influence the ground state realized in materials with 
various competing phases. Therefore, it will be useful 
to understand the competing phases and the transitions 
between them, based on symmetries that are broken by 
the small perturbations. 

In this work, we develop an effective Z 2 gauge the¬ 
ory, or equivalently a quantum dimer model (QDM), 
for competing phases in kagome antiferromagnets. The 
effective Z 2 gauge theory enables us to study compet¬ 
ing Z 2 SL and VBS phases on equal footing by consid¬ 
ering symmetry-allowed interactions in the correspond¬ 
ing QDM.®®’®® Among previous QDM approaches on 
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(a) 12-VBS (b) 36-VBS (c) 6-VBS 





FIG. 1. (Color online) Dimer configurations of the 12-site, 36-site, and 6-site VBS phases, (a) The 12-site VBS is characterized 
by the diamond-shaped dimer patterns (marked with light blue), (b) The 36-site VBS has dimer patterns of the perfect 
hexagons and stars (marked with light blue), (c) The 6-site VBS features parallel dimers (shaded with pink) and the zigzag 
dimers (shaded with green). 


kagome antiferromagnets,^®’^^’^®^^^ our theory is con¬ 
structed based on the recent Z 2 gauge theory by Wan and 
Tchernyshyov in Ref. [29]. They used the effective Z 2 
gauge theory to understand the topological Z 2 SL phase 
and its leading valence bond correlations, identified in the 
DMRG computations^^ on the nearest-neighbor Heisen¬ 
berg model on the ideal kagome lattice. Here we gen¬ 
eralize this theory to understand possible VBS phases 
and also extend the theory to incorporate the lowered 
rotational symmetry motivated by the vanadium oxyflu- 
oride compound. With this extension, we have another 
parameter for studying nearby VBS and Z 2 SL phases. 

On the ideal kagome lattice, we find that our Z 2 gauge 
theory captures not only the 12-site VBS but also the 
36-site VBS mentioned earlier [see Fig. 1 (a) and (b)]. 
In addition, we discover another interesting phase, the 
6-site VBS [Fig. 1 (c)] which is closely related to a par¬ 
ticular feature of the Z 2 SL identified in DMRG calcula¬ 
tions, i.e. symmetry-breaking valence bond modulations 
on certain cylinder geometries.These results show 
a remarkable consistency with the previous work consid¬ 
ering a different form of Z 2 gauge theory for the ideal 
kagome lattice.Almost all VBS patterns found in the 
latter are also identified in our Z 2 gauge theory. Gombin- 
ing our results on the VBS phases with the results in Ref. 
[29] on the Z 2 SL, we show that our gauge theory pro¬ 
vides a simple but generic model that captures various 
competing phases such as the Z 2 SL phase, the 12-site, 
36-site, and 6-site VBS phases. 

When the lowered symmetry by inequivalent up- and 
down-triangles is incorporated in the theory, we observe 
interesting effects of this on the VBS phases found on 
the ideal kagome lattice. First, the 12-site VBS is gen¬ 
erally stable to this perturbation. In the 36-site VBS, 
however, we find significant dimer melting, particularly 
at the “star” dimer structures while the perfect hexagon 


patterns remain relatively robust [the star and perfect 
hexagon dimer structures are marked with light blue in 
Fig. 1 (b)]. Interestingly, the 6-site VBS does not appear 
in the lowered symmetry environment. In the later part 
of this paper, we will discuss the implications of these 
findings and also compare the results with a different 
type of Z 2 gauge theory. 

A. Hamiltonian and overview 

Here we introduce our model Hamiltonian and provide 
an overview of this paper. For an intuitive picture, we 
start with a quantum dimer model on the kagome lattice 
that is equivalent to the Z 2 gauge theory. 

The QDM is defined in the Hilbert space of spin-singlet 
(or dimer) product states. Each state satisfies the so- 
called hardcore dimer constraint, i.e. each site on the 
lattice is covered by only one dimer. In the model, mo¬ 
tions and interactions of the dimers are described by the 
following Hamiltonian: 

Hqdm =—h Y.mD\+KY,eu\D){D\ 

+AY,Vt\T){T\. ( 1 ) 

The first term with the coupling h describes dimer mo¬ 
tions along various transition graphs (closed loops along 
which the dimers move) around each hexagon plaquette 
on the kagome lattice. The 32 transition graphs are listed 
in Table 1. The operator \D){D\ generates dimer motions 
between dimer configurations \D) and \D). As an exam¬ 
ple, Fig. 6 (a) shows the two dimer configurations on the 
12-length transition graph. The set ^ consists of all the 
dimer configurations included in the transition graphs. 
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Multiplicity 
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1 

3 
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3 

-3 

8 
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6 

-1 
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1 

10 

\ 

6 

-1 

10 

\“/ 

6 

1 

10 
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3 

3 

12 

V 

1 

-3 


TABLE I. (Color online) Transition graphs for the dimer mo¬ 
tions and values of the interaction energy coefficient eo- The 
transition graphs are denoted with light blue lines. In each 
graph, one of the dimer configurations (|D)) is depicted with 
red lines. The third column counts the number of equivalent 
graphs for a given shape. 


The second term with the coupling K represents dimer 
interactions with the energy coefficient of the con¬ 
figuration \D). Values of eu, listed in the last column 
of Table I, are determined from the Z 2 gauge theory 
as we will show later. Depending on the sign of K, 
the dimer interactions energetically favor different con¬ 
figurations; for K > 0, the diamond- and star-shaped 
dimer configurations are favored, while for K < 0, 
the hexagon- and doubled-diamond-shaped configura¬ 
tions are preferred [see Fig. 6 (b) and (c)]. The Hamil¬ 
tonian Hqdm with K > 0 and A = 0 corresponds to the 
QDM description of the Z 2 gauge theory considered by 
Wan and Tchernyshyov. 

The third term with the coupling A plays the role of a 
dimer potential energy. This dimer potential term incor¬ 
porates the lower symmetry as observed in the vanadium 
oxyfluoride compound by assigning different energy val¬ 
ues to the triangles with energy coefficient ryy as listed 
in Table II. The dimer potential energy rjT is determined 
from the Z 2 gauge theory shown later. If A > 0, dimers 
tend to occupy the up-triangles while leaving the down- 
triangles empty and vice versa. Therefore this potential 
energy term reflects the triangle inequivalence. 

Although it is not so obvious in this form, the QDM 
(and the equivalent Z 2 gauge theory) possesses two ex¬ 
actly solvable cases: {i) K = A = 0 and (ii) h = 0. In 
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TABLE II. (Color online) Values of the potential energy coef¬ 
ficient rjT for various dimer configurations on up- and down- 
triangles. Dimers are denoted with red lines. 


the former case, the ground state is exactly expressed by 
the RVB wave function on the kagome lattice, i.e. the 
superposition of all dimer configurations with equal am¬ 
plitudes, as pointed by Misguich, Serban, and Pasquier 
in Ref. [33]. This special point represents a Z 2 spin liq¬ 
uid phase. In the latter case, the ground state manifold 
contains various degenerate dimer configurations (singlet 
product states), each of which is an eigenstate of Hqum 
and corresponds to a valence bond solid phase. For in¬ 
stance, when K > 0 and h = A = 0, we find the singlet 
product state representing the 12-site VBS [Fig. 1 (a)] in 
the ground state manifold. When K < 0 and h = A = 0, 
the singlet product states corresponding to the 36-site 
VBS and 6-site VBS [Fig. 1 (b) and (c)] are found. 

In this paper, we explore the intermediate region be¬ 
tween the two exactly solvable cases. For convenience, we 
use the following parametrization for the coupling con- 


v-v>h 



FIG. 2. (Color online) Schematic phase diagram of the Z 2 
gauge theory. The exactly solvable points are denoted by 
the blue dot {K = A = 0; ■!? = 0) and the green line 
(/i = 0; = 7r/2). The red star represents the Z 2 SL ground 
state of the spin-1/2 nearest-neighbor Heisenberg model in 
the DMRG simulations, ip —>■ —ip takes A —>■ —A, which is 
equivalent to switching the up and down triangles. Without 
loss of generality we only study ip > 0. 
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stants. 

h = cos'd, 

K = sinz? cos(p, (2) 

A = sini? sin(/j, 

where 0 < '(9 < 7r/2 and —tt < ip < tt. With this 
parametrization, we draw a schematic phase diagram 
in Fig. 2. The exactly solvable cases are denoted by 
a blue dot [K = A = O; = 0) and a green line 
{h = 0; "d = 7r/2). The dashed line represents a schematic 
phase boundary separating the Z 2 SL and VBS phases 
near the exactly solvable points. In the phase diagram, 
the schematic location of the Z 2 SL ground state of the 
nearest-neighbor Heisenberg model is represented by a 
red star. In this study, we explore the VBS regime in the 
framework of the SL-to-VBS transition. 

In SL-to-VBS transitions, the central objects are the 
visons, topological vortex excitations in the Z 2 spin liquid 
phase.®’®® The transition is described by the condensa¬ 
tions of visons.The Z 2 gauge theory provides 
a convenient platform for the description of the visons. In 
Sec. II, we formulate the QDM on the kagome lattice by 
a Z 2 gauge theory on an effective honeycomb lattice. To 
describe the visons as local objects, the Z 2 gauge theory is 
transformed into the dual Ising model on a triangular lat¬ 
tice via a dual mapping in Sec. III. The transformations 
taken for our low energy theory are summarized in Fig. 
4. The SL-to-VBS transitions as well as the aforemen¬ 
tioned exactly solvable points will be investigated mainly 
based on the dual Ising model. 

In Sec. IV, we investigate VBS phases appearing from 
the transitions via Ginzburg-Landau functionals with vi- 
son condensation order parameters. The Landau theories 
are constructed based on projective symmetry group®® 
(PSG) analysis for the visons and soft mode expansion of 
the vison field. From the Landau theories and other com¬ 
plementary computations, we find 12-site VBS phases 
around the ~ 0 line and 36-site and 6-site VBS phases 
around the ~ tt line in the phase diagram (Fig. 2). 
Readers who are interested in the resulting VBS phases 
rather than the computations are advised to jump di¬ 
rectly to Sec. IV. 

In Sec. V, we compare our results with the previous 
study in Ref. [27] and discuss influence of the triangle 
inequivalence, incorporated by the dimer potential en¬ 
ergy A, on the 12-site, 36-site, and 6-site VBS phases. In 


vison condensation 


Z2 SL VBS 

FIG. 3. Transition from a Z 2 spin liquid to a valence bond 
solid by vison condensation. 
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FIG. 4. Transformations in the low energy effective theory. 

addition, we discuss the 6-site VBS as a promising low 
energy singlet fluctuation in the kagome antiferromag- 
net and make an interesting remark on it in relation to 
the recent DMRG simulation results. We provide several 
computational details in appendices. 

II. Z2 GAUGE THEORY 

Now we show the explicit construction of our Z 2 gauge 
theory. Equivalent to the quantum dimer model Intro¬ 
duced In the previous section, the gauge theory can be 
described by a simpler Hamiltonian in which the vison ex¬ 
citations are easily identihable. We transform the QDM 
Into a gauge theory by employing the dimer-arrow map¬ 
ping introduced by Elser and Zeng.®® 

Dimers on the kagome lattice are mapped to arrow rep¬ 
resentations defined on the honeycomb lattice obtained 
by connecting the centers of the kagome triangles [see 
Elg. 5 (a)]. The dimer-arrow mapping rules are depicted 
in Fig. 5 (b). The two-in-one-out arrow configuration 
corresponds to a dimer on the kagome link with the two 
incoming arrows (upper figure). On the other hand, the 
all-out arrow configuration means there is no dimer on 
a given kagome triangle (lower figure). By the mapping 
rules, we establish a one-to-one correspondence between 
the dimer configurations satisfying the hardcore dimer 
constraint and the arrow representations. 

This can also be expressed by the Z 2 variable cr“ (= 
±1), assigned to each link on the honeycomb lattice. If 
the arrow on a link ij goes from the sublattice A to B 
(B to A), we set ctT = -|-1 {afj = —1). With this Z 2 
representation, the hardcore dimer constraint takes the 
following form. 

^ _ X X X _ / Tl (f G 4I) 

Qi — — I 

where the sites, 1, 2, 3, represent the three nearest- 
neighbors of a given site i on the honeycomb lattice [Fig. 
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FIG. 5. (Color online) Dimer-arrow mapping for the Z 2 gauge 
theory formulation, (a) The honeycomb lattice, where the 
arrow representations are defined, is drawn in thick lines. The 
two sublattices, A and B, are denoted by solid and empty dots, 
(b) The arrow representations for a dimer-occupied triangle 
(upper) and an empty triangle (lower). A dimer is denoted 
by a red line. 


5 (b)]. 

Dynamics of the dimers can be incorporated by intro¬ 
ducing Z 2 variables conjugate to {crfj}. Both vari¬ 

ables have Pauli matrix representations as their notations 
imply. Kinetic terms encoded in the QDM [Eq. (1)] can 
be represented by a product of variables around a 
honeycomb plaquette a: 

0 ) 

{ij)ea 

This plaquette term encapsulates all the dimer motions 
described by the 32 transition graphs in Table I [for an 
example, see Fig. 6 (a)]. Notice that [F^,Qi] = 0 as the 
dimer motions respect the hardcore dimer constraint. 

In the Z 2 gauge theory description, the QDM takes the 
following form. 

Hzgt = -hY,F^^+K Y. ^E4- (5) 

“ {^3)//{kl) {ij) 


The first sum represents the dimer motions, and the sec¬ 
ond and third sums respectively correspond to the dimer 
interaction and potential energies in the QDM. In the sec¬ 
ond sum (K), each bilinear of is defined on parallel 
links, {ij)//{kl), of a honeycomb plaquette as depicted 
in Fig. 5 (a). It is straightforward to check that the sum 
of the bilinears is equivalent to the dimer interactions in 
the QDM [some cases in Table I are shown in Fig. 6 (b) 
and (c)]. It is also easy to check that the third sum (A) 
represents the dimer potential energy, e.g., from Fig. 5 
(b). We note that the Hamiltonian Hzgt with K > 0 
and A = 0 is the phenomenological model for the Z 2 spin 
liquid studied in Ref. [29] 

In this theory, we define the dimer occupation number 
d on the kagome lattice in the following way. 


di2 = + J = 


1 (occupied) 
0 (empty) 


( 6 ) 





FIG. 6. (Golor online) Descriptions of the h and K terms 
in the Hamiltonians, Eqs. (1) and (5). (a) One example of 
the dimer resonance motions (along the 12-length transition 
graph in Table I) by the h term, (b) The diamond- and star¬ 
shaped dimer patterns are favored by the positive K dimer 
interactions, (c) The hexagon- and doubled-diamond-shaped 
dimer patterns are favored by the negative K dimer interac¬ 
tions. 


Here, the dimer occupation number d\2 is defined on the 
kagome link that the honeycomb links il and i2 cross 
[see Fig. 5 (b)]. It can be easily checked that the above 
definition applies to both sublattices (t € A or H). Note 
that the average dimer occupation is d = 1/4 for any 
hardcore dimer configuration. 

Now we interpret the Hamiltonian Hzgt in tbe lan¬ 
guage of the Z 2 gauge theory. The Hamiltonian has the 
Z 2 gauge fields and conjugate electric fields {crfj} 

as underlying degrees of freedom. It is written in terms 
of two gauge-invariant operators, F^ and afj. These op¬ 
erators are invariant under the Z 2 gauge transformation: 
afj —>■ SiafjSj with Sij = ±1. The operator F* measures 
the Z 2 magnetic flux through a honeycomb plaquette, 
e.g., 0 flux for F^ = 1 and tt flux for F* = — 1. The elec¬ 
tric fields are generated by the Gauss law constraint, Eq. 
(3), with the alternating background Z 2 charges {Qi} on 
the honeycomb lattice. 

It is well known that the Z 2 gauge theory has two dis¬ 
tinct phases: the deconfined phase when h ^ |F|,|A| 
and the confined phase when h |F|, |A|.®^ In the de¬ 
confined phase, the ground state is characterized by a 
uniform 0 flux and its elementary excitations are the 
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gapped magnetic excitations carrying a tt flux called the 
visons.^^ As depicted in Fig. 7 (a), the visons can be 
created by the string operator 

n <’»• m 

{ij)es 

where the string s extends from the vison position to in¬ 
finity. The visons are obviously nonlocal objects in the 
gauge theory. As the electric terms {K and A) dominate 
the magnetic term {h), the vison energy spectra get more 
dispersive due to vison hoppings induced by the electric 
terms. When the vison gap is closed, vison condensa¬ 
tion leads to the confined phase accompanied by lattice 
symmetry breaking. In terms of the spin model, the de- 
confined and confined phases correspond to the Z 2 spin 
liquid and valence bond solid phases, respectively. 


III. DUAL ISING MODEL 


In this section, we map the Z 2 gauge theory to the 
dual Ising model in which the visons become local ob¬ 
jects. The dual Ising model enables us to identify the 
exactly solvable points of the theory and calculate the 
vison dispersions via a soft spin approximation. We con¬ 
sider four cases that represent different regions of the 
phase diagram (Fig. 2). For each case, we construct a 
Ginzburg-Landau functional of the vison fields in Sec. 
IV, from the projective symmetry group analysis of the 
visons. 

The Z 2 gauge theory on the honeycomb lattice is now 
transformed into the dual Ising model on a triangular 
lattice via the following dual mapping. 


pz _ X 
^ OL OL'! 

rrX _ \ 

^ij 


( 8 ) 


In this mapping, the u-operators on the honeycomb lat¬ 
tice are expressed in terms of the new r-operators on the 
dual triangular lattice [thick lines in Fig. 7 (b)]. The 
r-operators satisfy the same Pauli matrix algebra as the 
cr-operators. Notice that we are using Latin indices for 
the honeycomb lattice and Greek indices for the triangu¬ 
lar lattice. In this dual description, is the vison density 
operator and is the vison creation/annihilation oper¬ 
ator. Note that and the string operator in Eq. (7) are 
essentially the same up to a ± sign, and the visons are 
local entities on the dual triangular lattice. Dynamics of 
the visons are now governed by the dual Ising model: 


HdIM — ~h 'y [ T^-\-K 'y [ Xaj3\firyTy^T^-\-A y [ XoipT^T^ 
“ <((“.7))> 

( 9 ) 

Here, the K and A terms provide vison hoppings to the 
third and first nearest-neighbors on the triangular lattice, 
respectively. 

The additional variable A (= ±1) in the mapping is 
employed to keep the Gauss law constraint Eq. (3), which 





FIG. 7. (Color online) (a) Pictorial description of a vison. (b) 
The dual triangular lattice is depicted in thick lines and its 
sites are labeled with Greek indices, (c) The gauge choice for 
{Aa^} on the dual triangular lattice. The solid and dashed 
lines represent the links with Ac/s = -fl and A^/s = —1, re¬ 
spectively. The shaded rectangle indicates a unit cell in the 
gauge choice with two sublattices, 1 and 2. The blue arrows 
are the lattice vectors: u = (1, 0) and v = (0, %/3). 


now takes the following form. 


Q^ 




— 1 (a/Sy = A) 
-hi (a/37 = v) 


( 10 ) 


The constraint basically tells us the Berry phase that a 
vison will acquire when it hops around an elementary 
triangle. It must be noted that the dual mapping is in¬ 
variant under the Z 2 gauge transformation: 


Aa/3 GcG^Aq^, 


( 11 ) 


with Gq, = ±1. This fact allows many choices for the 
A-variables. We choose the gauge depicted in Fig. 7 (c), 
which is the same gauge choice taken in Ref. [29]. 


A. Exactly solvable points 

The dual Ising model has two exactly solvable cases: 
•(i) K = A = 0 in the deconfined regime and (ii) /i = 0 in 
the confined regime. In the former, the ground state is a 
paramagnetic state with = 1 at every site. In terms 
of the QDM, it corresponds to the superposition of all 
the hardcore dimer states with equal amplitudes, which 
is the resonating valence bond spin liquid state. In the 
latter case, the model has a degenerate ground state man¬ 
ifold with various hardcore dimer states, or valence bond 
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solid phases, depending on the angle ip [Eq. (2)]. When 
{'d,ip) = (7r/2,0) the degenerate ground states are all 
symmetry-related states of the 12-site VBS phase. There¬ 
fore the 12-site VBS is the only phase in the ground state 
manifold. In contrast, when ("d, ip) = (7r/2, tt) the ground 
state manifold contains various distinct VBS phases in 
addition to the 36-site VBS. The 6-site VBS in Fig. 1 (c) 
is one such phase. We find that the 12-site, 36-site, and 
6-site VBS phases appear in a finite region of the param¬ 
eter space ip-. 12-VBS in 0 < v? < ipa, and 36- and 6-VBS 
\n iph < ip < TT. We determine the phase boundaries, 
ipa = 1.222 and ipb = 1.571, by conducting a full search in 
a 24-site triangular cluster with periodic boundary con¬ 
ditions (which corresponds to a 72-site kagome cluster). 
Its results are summarized in Fig. 12, which shows the 
ground state energy and degeneracy as functions of ip. 
More details on the h = 0 case are provided in Appendix 
A. 

The phase diagram in Fig. 2 shows the exactly solvable 
points with the blue dot {K = A = 0) and the green line 
{h = 0). The intermediate regime between the two cases 
will be explored in Sec. IV. We shall find the 12-site, 
36-site, and 6-site VBS states as the broken symmetry 
phases that appear via vison condensation from the 'Ll 
SL (as we move outwards in the phase diagram in Fig. 2). 

B. Soft spin approximation: vison excitations 

Close to the transition point between SL and VBS, it is 
reasonable to take the following soft spin approximation: 

= il —y (j)a G (12a) 

. (12b) 

OL Ol 

Here, the Ising spin is replaced by a real-valued, 
coarse-grained vison field (pa , and effects of the transverse 
field term r® is described by the conjugate momentum 
i\<Pa,T^p\ = i^ap) and a vison mass m. The approxi¬ 
mation leads to the soft spin Hamiltonian: 

Hsoft = 2 ^ ^ ip^a ™ ^q) V 'y ) MapPa^Ppi (13) 
Oi 

where the matrix Map contains the hopping amplitudes 
from the K and A terms. Diagonalizing the Fourier 
transformed Hamiltonian, we obtain two bands of the 
vison dispersions, a;±(q), owing to the two sublattices. 

w±(q) = + Kfiq) ±\A\^/6 + /(q), (14) 

where /(q) = 2cos(2q • u) -|- 4sin(q • u)sin(q • v). 

Positions of the soft modes (the lowest energy modes) 
in q-space are marked with red dots in the figures in Ta¬ 
ble HI. Depending on the parameters, K (oc cosip) and 


case 1, 3 case 2, 4 



TABLE III. (Color online) Positions and degeneracies of the 
soft modes in the four cases. The soft mode positions are 
marked with red dots in the first Brillouin zone (rectangle) 
of the soft spin Hamiltonian Hsoft', Q = ~ 

^(1, %/3), and Q 2 = ^(5, %/3). The hexagon indicates the 
Brillouin zone of the kagome lattice. The green dots show 
positions of the lowest two-vison excitation computed by the 
equation: if2(q) = mink[ta-(k) -|- a;_(q — k)]. In the above 
table, the fourth column, Af, is the total number of the soft 
modes in each case. 


A (cx sin:^), the soft modes have different positions and 
degeneracies as summarized in the table. We consider 
four cases: (i) = 0, (ii) ip = tt, (hi) 0 < |i^| < and 

(iv) I < < TT. In cases 1 and 3, four soft modes occur 

at the points, ±Q, with a twofold degeneracy at each 
point (left figure). Case 2 has eight soft modes at the 
points, ±Qi^ 2 , each twofold degenerate, whereas case 4 
has four non-degenerate modes at the same points (right 
figure). In the remaining f ^ ly’l f, which is not in¬ 
cluded in the table, soft modes occur at incommensurate 
positions, which vary with p and have line degeneracies. 
Due to this complexity, the last case is not considered in 
our study. For each case in the table, Landau theory will 
be developed based on the vison soft modes. 

It is important to note that the visons are gauge- 
dependent objects owing to the gauge freedom in Eq. 
(11). Hence, the vison follows a projective representa¬ 
tion of the symmetry group (PSG) of the system and the 
soft mode positions can change depending on the gauge 
choice. When periodic boundary conditions are imposed 
on the system, the visons are always created in pairs by 
physical operators since Wa^a = 1- Hence the lowest 
two-vison excitations are useful in getting physical in¬ 
formation about the VBS phase stabilized by the vison 
condensation. From the lowest modes of the two-vison 
excitations (green dots in Table HI), we can figure out 
the size of the unit cell and the lattice structures of the 
VBS phases. Specifically, the VBS phases have triangu¬ 
lar lattice structures with 12-site unit cells in cases 1 and 
3, and also triangular lattice structures but with 36-site 
















unit cells in cases 2 and 4. Detailed VBS structures are 
studied via Landau theories later. 


C. Vison PSG 


The Af xAf matrix, Xrnn, is the PSG representation of X 
in the order parameter space {ipn}- Landau theories in 
the next section hinge on the above soft mode expansion 
and the vison PSG. 


Now we analyze the projective symmetry group of the 
visons in the four cases in Table III. Gases 1 and 2 have 
the space group symmetry SGi ^2 of the ideal kagome lat¬ 
tice. 

SGi,2 = Span{T„4,i?6}, (15) 

where the three generators represent the u-translation 
T„, the y-axis-reflection ly, and the 7r/3-rotation Rq (see 
Fig. 8). On the other hand, cases 3 and 4 have the 
symmetry group 

SG3.4 = Span{T„4,i?3}. (16) 

Here, the rotation symmetry is lowered to the 27r/3- 
rotation i ?3 due to the triangle-inequivalence by a 
nonzero A term. The visons (((>q) transform projectively 
under the symmetry operation. The projective sym¬ 
metry group,®® PSG = {GxX}, is an extended sym¬ 
metry group with the Z 2 gauge transformation Gx for 
each symmetry operation X (here, the invariant gauge 
group is IGG=Z 2 ). The gauge transformation Gx is 
determined by the condition that the vison Hamilto¬ 
nian is invariant under the PSG, i.e. Mx(a)x(p) = 
AIapGx[X{a)]Gx[X{l3)]. We describe the generators of 
the PSG pictorially in Fig. 8 and put their explicit ex¬ 
pressions in Appendix B. 

Next, we turn our attention to the vison condensation 
order parameters and discuss how to extract a representa¬ 
tion of the PSG in the order parameter space. The order 
parameters are defined by expanding the vison field in 
soft modes as the following. 

= (17) 

n—1 

Here, r and s (= 1,2) are the unit cell and sublattice 
indices [Fig. 7 (c)], and N is the number of soft modes 
available. The complex number, ipn, is the vison order 
parameter corresponding to the n-th soft mode with wave 
vector q„ and eigenvector v” of Hgoft- With the above 
soft mode expansion, we deduce the PSG of the order 
parameters as follows. 

GxX[Mr)] = 

n 

= (18) 
n 

In the first line, (r',s') = X(r, s). The PSG is read 
off from the second equality in the form of the matrix 
equation: 

'4’m ~ ^ ^ Xran'<Pn- (19) 


IV. LANDAU THEORIES FOR VALENCE 
BOND SOLIDS 

We now construct Landau functionals for the four cases 
in Table HI. Recall that the four cases represent differ¬ 
ent regions with different space group symmetries (SGi ,2 
and SG 3 _ 4 ) in the phase diagram in Fig. 2. The Lan¬ 
dau theories enable us to investigate possible symmetry 
breaking patterns that can appear in the valence bond 
solid phases of the dual Ising model. We discuss the 
symmetry breaking patterns at mean-field theory level. 


A. Case 1 : = 0 

Case 1 has four vison soft modes at ±Q with the sym¬ 
metry group SGi_ 2 - Accordingly, the vison field expan¬ 
sion is 


0(r) = -I- -I- c.c. (20) 


with eigenvectors, = (1,0)^ and = (0,1)^, and 
order parameters, = (^/>i,^ 2 ,V’l,f/’^)- The PSG rep¬ 
resentation carried by the order parameters is generated 
by the three matrices: 


% = 


r — 
-^v ~ 


i 

0 

0 

0 

0 

i 

0 

0 

0 

0 

—i 

0 

0 

0 

0 

—i 

0 

0 

1 

0 \ 

0 

0 

0 

-M 

1 

0 

0 

0 

0 

i 

0 

0 / 



/ 1 -f 

1 -i 

i 

-1 -i 

1 i 

2 

1 i 

1 i 


\ 1 -i 

— 1 i 


(21a) 


(21b) 


(21c) 


These generate a 48-element subgroup of 0(4) isomor¬ 
phic to GL(2,Z3)®®. With the PSG representation, we 
construct a fourth order Landau functional that is in¬ 
variant under this group. 


£ = |aT|2 + + aih + 02 X 2 (22) 

with 

Xi = —&p\p\ + p^cos4di -I- p2Cos 402, (23a) 

X 2 = /OiP2sin20i [cos(0i -I- 62 ) + sin(di -|- 62 )] 

+ pip2sin202[cos(0i -I- 62 ) - sin(di -f 62 )] 

-I- p\p\s\n[2{9i - 02)]. (23b) 
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FIG. 8. (Color online) Pictorial descriptions of the generators in the vison PSG. (a) The generators of the symmetry group 
SG 3 , 4 . Tu- translation by the vector u. ly-. reflection with respect to the y-axis. iZs: rotation by 27r/3 with respect to the 
origin. For the symmetry group SGi, 2 , the f? 3 -rotation is replaced with the rotation by tt/S (Re)- (b,c,d) The action of the 
symmetry operation, X{g {ly, Ra, R3}), on the Hamiltonian Hsoft with the gauge choice {Aa/s} in Fig. 7. Red dots in each 
figure indicate the sites where the sign changes in the associated gauge transformation Gx ■ The translation operator Tu is not 
shown here due to its trivial action. 


Here we have used the parametrization ipn = Pne*®"(n = 
1,2) with real-valued pn and 0„. We will be using the 
same parametrization for the other cases as well. 

With the Landau functional, symmetry breaking pat¬ 
terns in VBS phases are investigated via a mean-field 
approach in the confined regime (ik 7 ^ 0). We distinguish 
different phases by observing changes in the mean-field 
energy profile and the dimerization pattern. The dimer¬ 
ization pattern is identified by computing the dimer oc¬ 
cupation number d [Eq. ( 6 )] at each kagome link; d is cal¬ 
culated by using the dual mapping rule [Eq. ( 8 )] and the 
soft mode expansion [Eq. (20) in this case] with mean- 
field solution 'I’ mf- In plotting the dimerization pattern. 



FIG. 9. (Golor online) Mean-field phase diagram of the Lan¬ 
dau theory [Eq. (22)] for case 1. VBSia phase is the 12-site 
VBS in Fig. 1 (a). The phase diagram is obtained with fixed 
u = 10 . 


we use red {d > d), black {d < d), and dashed {d = d) 
lines. Remember that c? = 1/4 is the average dimer occu¬ 
pation. In our plotting scheme, the thickness of the red 
and black lines is proportional to the deviation from the 
average, [d — J). 

Figure 9 shows a mean-field phase diagram in the 01-02 
plane of the above Landau theory. The diagram contains 
three different phases: VBSia, VBSib, and VBSic. 

• Phase A . The VBSia phase is the 12-site VBS fea¬ 
tured with the dimer structures of the diamonds 
and the stars [Fig. 1 (a)]. In this phase, the red 
dimers with d > d respect the hardcore dimer con¬ 
straint (each site is occupied by only one red dimer) 
since it is dominated by a particular singlet product 
state. 

• Phases B and C . In contrast to phase A, VBSis 
and VBSic phases consist of the black and dashed 
bonds without any red dimers. The absence of the 
red dimers implies that the phases are superpo¬ 
sitions of various dimer product states without a 
clear symmetry breaking pattern. Particularly, the 
VBSic has uniform dimer density d = d despite 
having a nonzero vison condensation mf 7 ^ 0 .®^ 
In this phase, symmetry breaking is anticipated to 
occur when fluctuation effects beyond the mean- 
held are incorporated. The VBSib may be under¬ 
stood as a variation from the VBSia with the star 
dimer structures being melted by dimer resonance. 


B. Case 2-. ip = -n 

In case 2 , there are eight soft modes at ±Qi .2 with the 
symmetry group SGi_ 2 - With the soft modes, vison held 
is expressed as 

<f>{v) = -b -b c.c. 

+ -b c.c. (24) 
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FIG. 10. (Color online) Mean-field phase diagrams of the Landau theory [Eq. (26)] for case 2. The phases VBS 2 A, VBS 2 a'i 
VBS 2 A", and VBS 2 A'" overall have the dimer structure of the 36-site VBS in Fig. 1 (b) with modifications in the dimer 
orientations in the hexagons and the overall dimer density. The phases VBS 2 fl, VBS 2 B', and VBS 2 B" have their origin in the 
6-site VBS in Fig. 1 (c) with different patterns of dimer density modulation. In each phase diagram, we set u = 10 and the 
other a-parameters to be zero ( 03,4 = 0 in the left and ai ^2 = 0 in the right). 


with the same eigenvectors as in the case 1 and or¬ 
der parameters 'll = ('*/’i,'03,"04; V’l;' 02 ;V’D- The 
order parameters have the following matrix representa¬ 
tion of the vison PSG. 


where z = These matrices generate a 288-element 

subgroup of 0(8) isomorphic to GL(2, Z 3 ) xD(3)®^. With 
the PSG representation, the following fourth order Lan¬ 
dau functional is obtained. 


Tu = 


r — 
-^y ~ 


TZe = 
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0 0 
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0 

0 
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0 0 
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z4 
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0 
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0 
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0 

0 

1 

z ® 

1 

z® 

0 
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1 

0 

0 z® 

z " 

z 2 

z 

0 
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(25a) 


(25b) 


,(25c) 


£ = | aT -|2 + r |«'|2 + u | T '| 4 +^ ( 26 ) 

m—1 

where 

T 22, 22, 22, 22 

Ti — P1P2 + P2P3 + P3P4 + P4P1 

- p\pI [1 -b cos2(di -b 6»3)] 


- P2P4 [1 + cos2(d2 + 6*4)], ( 27 a) 

n- 2 2,22 

T2 — P1P3 + P2P4 

- PiP2/33/04COs(6>i - 6*2 - 6*3 -b 6*4) 

- piP 2 P 3 PiVisin{ 9 i - 6*2 - 6*3 -b 6*4), ( 27 b) 

T 3 = Ji + J 2 - 4(^3 - 1) J 4 , (27c) 

I 4 = -(^3 + 1)^2 + J 3 + 4 J 4 . (27d) 


The explicit forms of Ji^ 2 , 3 ,a are relegated to Appendix 
G due to their lengthy expressions. 

Mean-field phase diagrams in this Landau theory are 
shown in Fig. 10: in the 04-02 plane (left) and 03-04 
plane (right). The phases appearing in the diagrams are 
classified into three groups: A, B, and G. The VBS phases 
are named after the groups the phases belong to; VBS 2 A, 
VBS 2 A', VBS 2 A", and VBS 2 A'" are of type A. 

• Group A. The phases in this group represent the 
36-site VBS phases characterized by the dimer 
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structures of the stars and hexagons [Fig. 1 (b)]. 
Among the phases, detailed dimer structures vary 
in aspects of the dimer orientations in the hexagons 
and the overall dimer density. To be specific, in 
VBS 2 yi and VBS 2 A' two different dimer orienta¬ 
tions are found in the hexagons whereas only one 
orientation is found in the hexagons of VBS 2 A"- In 
terms of the dimer density, VBS 2 A' and VBS 2 A" 
have notable modulations in the dimer density 
across the red dimers while VBS 2 A has more or less 
uniform density in its red dimers. In the VBS 2 A'" 
phase, only the bridge dimers, the dimers con¬ 
necting the stars and hexagons, are observed and 
the other star and hexagon dimer structures have 
melted away. It is interesting to note that VBS 2 A'" 
is similar to VBSis in terms of the dimer melting. 

• Group B. The three phases of this group have their 
origin in the 6 -site VBS pattern composed of par¬ 
allel and zigzag dimers [Fig. 1 (c)]. Here, the unit 
cell is enlarged by a dimer density modulation to 
36 sites. The phases VBS 2 B, VBS 2 B', and VBS 2 B" 
have different modulation patterns. In VBS 2 B' and 
VBS 2 _b", the dimer density is concentrated on the 
parallel and zigzag dimers, respectively. In VBS 2 B, 
the modulation occurs in both types of dimers. 

• Group G. VBS 2 C has no clear dimerization pattern 
since the dimer density is d = d at every link. This 
is essentially the same phase as VBSic at the mean- 
field level. 


C. Case 3: 0 < y; < ^ 

Now we consider case 3 that has four soft modes at ±Q 
with the symmetry group SGsa- In this case, the soft 
mode expansion of the vison field is exactly the same as 
in case I given by Eq. (20). The only difference between 
the two cases is the symmetry group, which is now SGa^ 
with the lowered rotation symmetry R 3 . We find that 
the Landau functional with the symmetry group remains 
the same as that of case 1 up to sixth order. Therefore, 
the mean-field phases in case 3 are the same as those in 
case 1 to this order. This fact suggests that the 12-site 
VBS phase is robust against the triangle inequivalence. 


D. Case 4: f < < tt 

In case 4, there are four soft modes at ±Qi ,2 under 
the symmetry group SGs^. Landau theory for this case 
is developed by the following soft mode expansion: 

</>(r) = -b -b c.c. (28) 



fli 


FIG. 11. (Color online) Mean-field phase diagram of the Lan¬ 
dau theory [Eq. (31)] for case 4. In the four phases of the dia¬ 
gram, dimer melting effects, indicated by dashed lines {d = d), 
are observed in common. The VBS 4 A and VBS 4 S are consid¬ 
ered as variants of the 36-site VBS by the dimer melting. 
The same phase diagram is obtained in two different cases: 
(i) u = 0, V = 20, and (ii) r = 0, v = 10. 


Here the eigenvectors, take different forms from the 

previous cases: 


v(i) ^i_/^2-y3-bt\ ^(2) ^J_/2-by3-b*\ 

vi \ 2 y ’ 1^2 \ 2 y ’ 

(29) 

with the normalization factors, vi = 2-\/3 — -s/S and 1^2 = 
2 Vs -b -/S. The order parameters, T' = (V'l, ^/’ 2 , , ^” 2 ) j 
have the following PSG representation. 



0 

0 

0 ^ 



0 

z® 

0 

0 



0 

0 

z-^ 

0 

5 


VO 

0 

0 

Z-® y 



f 0 

0 

z^ 0 \ 


0 

0 

0 zii 



^7 

0 

0 0 

5 

V 0 


0 0 




(z 

5/2 

0 


0 z-^ \ 

1 


0 2 

- 11/2 

z2 0 

V2 


0 

z-4 

z 

5/2 0 


V ^ 

-2 

0 


0 zii /2 y 


(30a) 
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where z = The group generated is isomorphic to 

(C 3 X SL(2, Z 3 )) X C 2 ,®^ which is a finite subgroup of 0(4) 
with 144 elements. The PSG representation leads to the 
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following sixth order Landau functional: 

£ = | a ^|2 + + y |^|4 + ^|^|6 ^ ^ 

with 

2^1 = {Pi - pl)PiP 2 ^os 2 {ei + 6 * 2 ), (32a) 

I 2 = p\p 2 Cos{Wi - 62 ) + pip\cos{9i - 56»2). (32b) 

Figure 11 shows four different mean-field phases in the 
01-02 plane. In all of the four phases, dimer melting 
effects, indicated by dashed lines {d = d), are observed. 

• Phases A and B . ¥ 684 ^ and VBS 4 B are consid¬ 
ered as variants of the 36-site VBS. VBS 4 A pos¬ 
sesses hexagon and bridge dimer structures al¬ 
though the stars are all erased out. In VBS 4 B, most 
of the 36-site VBS structures are wiped out leaving 
a triangular arrangement of the hexagon structures. 

• Phases C and D . VBS 4 C and VBS 4 D have similar 
structures to those of VBS 2 A'" and VBSis. Red 
dimers are found scattered in the VBS 4 c. Interest¬ 
ingly, the red dimers do not fit to any characteristic 
dimer structure observed in the 12-site, 36-site, and 
6 -site VBS phases. 

V. DISCUSSIONS 

In this paper, we studied the Z 2 gauge theory repre¬ 
sented by the Hamiltonian Hzgt [Eq. (5)] as an effective 
theory for low energy spin-singlet fluctuations in kagome 
lattice antiferromagnets. Our study was focused on va¬ 
lence bond solid phases in the intermediate regime sand¬ 
wiched between the two exactly solvable points in the 
phase diagram of Fig. 2. It is found that the 12-site VBS 
phase with the diamond pattern is the only VBS phase 
that can be realized on the p = 0 line. This is supported 
by three different approaches: the ground state manifold 
of the Hamiltonian when h = 0, the wave vectors of the 
lowest two-vison excitations, and the dimerization pat¬ 
tern from the Landau theory (VBSia). According to re¬ 
cent DMRG simulations,^® the 12-site VBS is considered 
the long-range ordered counterpart of the Z 2 spin liquid 
ground state of the spin-1/2 nearest-neighbor Heisenberg 
model. 

On the other hand, two different types of VBS phases 
are found on the p = tt line of the phase diagram. One 
of the phases is the 36-site VBS that was proposed as 
a ground state candidate of the Heisenberg model based 
on various studies, including the quantum dimer model^® 
and the series expansion.We have identified three 36- 
site VBS patterns, VBS 2 A, VBS 2 A', and VBS 2 A"- The 
remaining phase is the 6 -site VBS. In our Landau the¬ 
ory, the phase appears with an enlarged 36-site unit cell 
by dimer density modulations as in the phases VBS 2 B, 
VBS 2 B', and VBS 2 S''. It is interesting to note that the 
6 -site VBS does not have at all the hexagon and doubled- 
diamond dimer structures favored by the dimer interac¬ 
tion when K < 0 [compare Fig. 1 (c) and Fig. 6 (c)]. It 


Our theory 

Ref. [27] 

VBSia 

VBS If 

VBSifl 

VBS 2f 

VBS2A,2A' 

VBS 1a 

VBS 2 A'" 

VBS 2a 

VBS2fl' 

VBS 4a 

VBS2B" 

VBS 3a 


TABLE IV. Correspondence between the results of this paper 
and Ref. [27] 


contrasts with the 36-site VBS that contains those two 
dimer structures. 

The above results are remarkably consistent with a pre¬ 
vious work.^^ In a different form of the Z 2 gauge the¬ 
ory, they proposed 12-site VBS and 36-site VBS states 
as promising VBS phases in confinement transitions from 
the Z 2 spin liquid phase. While the representations of the 
visons differ in the two theories, the group under which 
the visons transform are isomorphic to one another in 
each case. The ip = 0 {(p = ir) corresponds to the fer¬ 
romagnetic (antiferromagnetic) next nearest neighbor in¬ 
teraction considered in Ref. [27]. Also, as summarized 
in Table IV, there is a surprising correspondence between 
the VBS patterns. Interestingly, the 6 -site VBS patterns, 
VBS 2 B' and VBS 2 B", in our theory also appear in the 
other theory. 

The appearance of the 6 -site VBS state in the Landau 
theory approach as well as in the ground state manifold 
of Hzgt with K < 0 and h = 0 suggests that it may be 
a competing low energy state considered in the kagome 
antiferromagnets. In fact, it was hinted in the DMRG 
simulations of the Heisenberg model. On the families of 
cylinders, YG4m-|-2 and YG(4m-|-l) —2, the ground state 
of the Heisenberg model displays symmetry-breaking va¬ 
lence bond modulation patterns; strong valence bonds 
formed along zigzag lines and straight lines. The mod¬ 
ulation patterns are understood as a resonating state of 
the 6 -site VBS dimer configurations with the zigzag and 
straight lines as transition graphs [see Fig. 1 (c)]. 

Beyond the ideal kagome lattice geometry, we also in¬ 
vestigated effects of the triangle inequivalence on VBS 
order formation. Upon introducing inequivalent trian¬ 
gles via the dimer potential energy term A, the 12-site 
VBS and 36-site VBS respond in different ways. Around 
the line :/? = 0 of Fig. 2, the 12-site VBS remains stable 
to anisotropy as suggested by the unaffected (i) ground 
state manifold of the Hamiltonian with /i = 0, (ii) vi- 
son soft modes, and (iii) Landau theory with the lowered 
rotation symmetry (the case 3). 

Around the p = n line, the inequivalent triangles 
have dimer melting effects on the 36-site VBS phase as 
represented by the dimerization patterns in VBS 4 A and 
VBS 4 B. In both patterns, the star dimer structures are 
all melted whereas the hexagon structures remain at least 
partly. This trend can be understood by noting that the 
hexagon structures carry empty triangles energetically fa- 
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VBS 

size of unit cell 

Rs 

As 


lA 

12 

/ 

/ 


IB 

12 

/ 

/ 

/ 

1C 

3 

/ 

/ 

/ 

2A,2A' ,2A" ,2A'" 

36 


/ 


2B,2B',2B" 

36 




2C 

3 

/ 

/ 

/ 

AA,AB,AC, AD 

36 


/ 

/ 


TABLE V. Remaining symmetries in the VBS phases. 


vored by the dimer potential A while the star structures 
do not have any empty triangles. Hence, the hexagons 
are more stable. At the stars, dimer resonances, which 
lead to dimer melting, are enhanced by the inequivalent 
triangle geometry. Interestingly, this is in contrast to the 
stability of the 12-site VBS (triangular array of the stars) 
around the (p = 0 line. The different behaviors at ~ 0 
and ~ TT are attributed to their different structures of 
low energy space of the Hamiltonian. We also note that 
no analogue of the 6-site VBS is found in the inequivalent 
triangle geometry. 

We have also attempted to investigate effects of the tri¬ 
angle anisotropy by extending the previous Z 2 gauge the¬ 
ory in Ref. [27] . As discussed earlier, both versions of Z 2 
gauge theory capture the 12-site and 36-site VBS phases 
on the ideal kagome lattice. As we tune anisotropy, both 
theories predict regions of 12-site VBS, incommensurate 
VBS, and 36-site VBS phases. They both show similar 
symmetry breaking patterns by the triangle anisotropy 
on the 36-site VBS phase. However, the 12-site VBS is 
destroyed immediately in favor of the incommensurate 
phase in the conventional theory, while it survives up to 
a finite critical anisotropy in our theory. The latter dif¬ 
ference is ascribed to difference in the way the triangle 
anisotropy is incorporated in each theory. Details on the 
conventional Z 2 gauge theory with triangle anisotropy 
can be found in Appendix D. 

Now we discuss remaining symmetries in the VBS 
phases found in this work. For each phase, size of unit 
cell and remaining point group symmetries are summa¬ 
rized in Table V. In the 12-site VBS phase (lA), the 
i?6 rotation symmetries with respect to the centers of 
the dimer stars are preserved while the ly inversions are 
completely broken. In comparison with the 12-site VBS, 
the 36-site VBS phases (2A, 2A', 2A") have less symme¬ 
tries due to more complicated and larger unit cell: only 
the i ?3 rotations with respect to the centers of the dimer 
stars among the point group symmetries. In other 36-site 
VBS phases (2H, 2B', 2H"), when considered as dimer- 
modulated 6-site VBS states, all the point group sym¬ 
metries are completely broken. In contrast to the above 
phases, the phases 1C and 2C may seem to have no bro¬ 
ken symmetries with uniform dimer density. However, 
these uniform dimer density phases are results obtained 
at the mean-held level. We expect symmetry breaking 
pattern would be emergent when huctuation effects are 


incorporated in those phases. In the anisotropic kagome 
lattice case, the VBS phases, 4A, AB, 4(7, AD, have the 
ly inversions as well as the i ?3 rotations. It is interesting 
to note that the VBS phases on the anisotropic kagome 
lattice preserve more symmetries compared to the 36-site 
VBS phases found on the isotropic kagome lattice. 

In this work, we considered a simple and generic quan¬ 
tum dimer model that could describe low energy singlet 
physics of kagome antiferromagnets. In regard to anti- 
ferromagnets on the anisotropic kagome lattice, it would 
be an interesting problem to investigate which micro¬ 
scopic spin models can be mapped onto our effective 
quantum dimer model. In the present situation with¬ 
out any numerical studies on anisotropic kagome lattice 
models, Mila’s effective model construction^^ on the J- 
J' anisotropic kagome Heisenberg model can be helpful 
for understanding the connection between the quantum 
dimer model and J-J' Heisenberg model. It would also 
be interesting to compare Mila’s effective model with our 
quantum dimer model on the anisotropic kagome lattice. 
We leave the investigation in this direction as a future 
study. 


ACKNOWLEDGMENTS 

We thank Robert Schaffer for helpful discussions and 
collaborations on a related study. This work was sup¬ 
ported by the NSERC of Canada, the Canadian Institute 
for Advanced Research, and the Center for Quantum Ma¬ 
terials at the University of Toronto. Some parts of this 
work (YBK) were performed at the KITP, supported by 
NSF PHYlI-25915. 


Appendix A: Exactly solvable points at h = 0 

Here we discuss in more detail the h = 0 (i? = 7r/2) 
exactly solvable points of the phase diagram in Fig. 2. 
At these points, the dual Ising model Hdim is readily 
diagonalized as follows. 

K E + A (Al) 

(((a.7»> (a./ 3 > 

with K = cosif and A = sintp. To explore the ground 
state manifold of the Hamiltonian, we put the system on 
a Hnite size cluster of the dual triangular lattice with a 
periodic boundary condition. 

We design the cluster to capture the 12-site and 36-site 
VBS states in Fig. 1. It is easy to see that the 12-site 
VBS state is the exact ground state at A = 0 by noting 
that the state consists of the star and diamond dimer 
configurations with the minimum energy of the K poten¬ 
tial [see Fig. 1 (a) and Table Ij. In the dual Ising model 
description, r'’’-spin configurations representing the 12- 
site VBS have a 8-site unit cell on the triangular lattice, 
which constrains the cluster to have a multiple of 8 sites. 









14 




FIG. 12. (Color online) Plots of the ground state energy of the 
Hamiltonian at /i = 0 and degeneracy as functions of the angle 
parameter tp. The results are obtained with a full search in the 
24-site cluster (at the top) with periodic boundary conditions. 
For explanations of the plots, see Appendix A. 


site is —K = —cost/? (green line in the figure). One can 
easily check the energy by acting the QDM Hamiltonian 
Hqdm with h = 0 on the dimer product state represent¬ 
ing the 12-site VBS [Fig. 1 (a)]. 


When </3h < < TT, the ground state manifold consists 

of various VBS phases including the 36-site VBS and 
6-site VBS. Depending on the presence/absence of the 
dimer potential energy (A) in the Hamiltonian, empty 
triangles in the VBS states are arranged in different ways. 
When iph < ip < TT, the empty triangles are all down- 
triangles in both 36-site and 6-site VBSs. The energy 
of these states, {K — A)/3 = {cosp — sin(^)/3, is plotted 
in blue in the figure. However, at the point </? = tt, the 
up- and down-triangles have no distinction due to the ab¬ 
sence of dimer potential. This fact leads to more choices 
on the distribution of the empty triangles in the VBS, 
which is reflected in the increased ground state degener¬ 
acy at that point (see the lower plot in Fig. 12). The 
pink line represents the energy of the 36-site VBS and 
modified 6-site VBS states with the same number of up- 
and down-empty-triangles: K/3 = cos:/9/3. The 36-site 
VBS is shown in Fig. 1 (b). The modified 6-site VBS 
is obtained by arranging the parallel dimers of the 6-site 
VBS [Fig. 1 (c)] in such a way that lines of the paral¬ 
lel dimers have alternating dimer configurations. Note 
that the empty-triangles are distributed along the paral¬ 
lel dimer lines in the 6-site VBS. 


In the intermediate region pa < < ‘fib, the dimer 

potential energy dominates the interaction energy, and 
none of the 12-site, 36-site, and 6-site VBS states are 
found. This region corresponds to the parameter region 
f ^ I'pI ^ f of Hsoft where vison soft modes occur at 
incommensurate momentum positions (Sec. HIB). 


On the other hand, possibility of the 36-site VBS state 
(strictly speaking, the size of unit cell) as the ground 
state was hinted in the soft-spin approximation approach 
(Sec. HIB). VBS states with a 36-site unit cell require 
the cluster to have a multiple of 12 sites. Combining 
the above two conditions, we find that the cluster size 
should be 24n with n = l,2,'--. As a cluster satisfy¬ 
ing this condition, we consider the 24-site cluster in Fig. 
12, which corresponds to a 72-site cluster on the kagome 
lattice. By performing a full search on the cluster, we 
obtain the ground state energy and degeneracy as shown 
in the figure (red dots as functions of the angle parame¬ 
ter (/?). In the parameter range 0 < :/? < tt, we find three 
different regions separated by two points, pa = 1.222 and 
Pb = 1.571. 

First, in the region 0 < </? < pa, the ground state mani¬ 
fold only contains dimer configurations of the 12-site VBS 
phase. In this case, the ground state energy per kagome 


Appendix B: Vison PSG 


In this appendix, we provide explicit expressions for 
the generators of the vison PSG in Fig. 8. First, physical 
symmetry part of each generator on the dual triangular 
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lattice is expressed as follows. 

T„ : (m, n)s (to + l,n)s, (Bla) 


f (TO,n)i -)> (-TO,n)i 

( (to, n)2 —>■ (—TO — 1, 71)2 ’ 


(Bib) 


i?6 : < 


(TO,n)i -)■ ( 

(to, 77)2 -)> ( 


m—3n m+n \ 

2 ’ 2 11 

m—3n—2 m+n 
2 ’ 2 


)2 


(TO,n)i ^ 

(TO,n)2^(«^,^^^)l 


(m + n : even) 

(to + n : even) 

? 

(to + n : odd) 

(to + n : odd) 

(Blc) 


R?. '■ < 


(TO,n)i 

(to, 72)2 

(to, 77)1 

(777,77)2 


f—7n—3n m—n\ 

\ 2 ’ 2 ii 

/ —m—3n—2 m—n\ 

\ 2 ) 2 71 

/—m—3n—1 m—n—1\ 

i— j —> — 2 — 12 

/ —m—3n—3 m—n—1\ 

\ 2 > 2 72 


(to + 77 : even) 
(to + 77 : even) 

(to + 77 : odd) 
(to + 77 : odd) 

(Bid) 


Here, the notation (to, 77 );; represents the sublattice site 
s (= 1 , 2 ) in the unit cell located at tou + ttv (to, 77 G Z). 
For the above symmetry operations, the accompanied 
gauge transformations can be chosen as follows. 


GT^(m,n)s 

Gijm,n)s 


Gue {m,n)s 


Grs {m,n)s 


n-ir (s = l) 
I (-!)”+! (s = 2) 
/(-ir (s = l) 

l(-l)” (5 = 2)’ 


(B2a) 

(B2b) 

(B2c) 



/ s = 1; mod(TO + 77, 4 ) = 0,1 \ 
s = 2; mod(TO + 77, 4 ) = 2, 3 J 
(otherwise) 


(B2d) 


Appendix C: Invariants I 34 in the case 2 

The invariants I 34 in Eq. (27 c,d) are composed of 
four terms Ji.2,3,4: 
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Ji = PiPsPa 


P 2 P 3 P 4 


P1P2P3 


cos(2^i + ^3 + ^4) + (—2 + '\/ 3 )sin( 2 ^i + ^3 H- ^4) 

—cos( 2^2 H“ ^3 “t“ O4) + (—2 + '\/ 3 )sin( 2^2 H“ ^3 “h O4) 

(—2 + '\/3)cos(^i -\- 02 2 ^ 3 ) + siii(^i 62 2 ^ 3 ) 


PiP2p\ (2 — '\/ 3 )cos(^i + ^2 H“ 2^4) + sin(^i -\- 62 2^4) 


J2 = {pIp2 - PipI) (2 - y 3 )cos( 6 »i - 62) + sin( 0 i - 62) 
- {pIpa - PspI) cos(6»3 - 0a) + (2 - y3)sin(6»3 - 04) 


J3 = pIpI r-v^cos 2 ( 6 li - 02 ) + sin 2 ( 6 li - 6 » 2 ) 

+ P3P4 v^cos2(03 - 0a) + sin 2(03 - 0a) 

+ iPiP 2 pl - PiP 2 pI) (-1 + y3)cos(6>i - 02 ) + (1 + '/3)sin(6»i - 02 ) 

+ (-PiP3P4 + P 2 /O 3 P 4 ) (1 + \/3)cos(6>3 - 6»4) + (-1 + V^)sin (03 - 6»4) , 

Ja = / 9 l/ 92 / 53 P 4 Sin( 6 »i - 02 + 03 - Oa)- 


(Cla) 


(Clb) 


(Clc) 

(Cld) 


Appendix D: Comparison with the conventional 
gauge theory 

Here we compare the results of the anisotropic kagome 
lattice with traditional Z2 gauge theory used in Ref. [ 27 ]. 
In this, the visons live on a fully frustrated dice lattice, 
corresponding to centers of the triangles and hexagons of 
the kagome lattice. (Fig. 13 ) As noted in Ref. [ 29 ], the vi¬ 
sons at the centers of triangles are higher in energy than 
the ones in the hexagons and the lowest energy descrip¬ 
tion can be done using the visons on a triangle lattice. 
The low energy results of both theories should be equiv¬ 
alent and is confirmed below in our case. 

The Hamiltonian of the visons can be written as 

n =, (Dl) 

where the couplings J^ ’s are chosen such that the prod¬ 
ucts around a plaquette are negative. With anisotropy, 
we have two independent nearest neighbor vison inter¬ 
actions, denoted Ji and J2. These correspond to in¬ 
teraction between visons living in the center of hexagon 
and up/down triangles. When Ji = J2 the vison bands 
are flat in momentum space, and in Ref. [ 27 ] ferromag¬ 
netic or antiferromagnetic next nearest neighbor interac¬ 
tions connecting two 3 -coordinated sites (t) were added 
to give dispersion to the flat vison bands. Ji 7^ J2 im¬ 
mediately gives dispersion to the visons and is smoothly 
connected to the antiferromagnetic t-interaction consid¬ 
ered in Ref. [ 27 ]. The states resulting from ferromagnetic 
next nearest neighbor interactions were immediately de¬ 
stroyed (which is different from the results of the main 
text) and the soft modes were found to be incommensu¬ 


rate to the lattice. This corresponds to states found in 
the main text at 7r/3 < < 7r/2. Here we only consider 

the t < 0 case. 


FIG. 13. (Color online) (left) Unit cell of a fully frustrated 
dice lattice with nearest neighbor interactions. It has 12 sites 
in the unit cell and the black (red) lines denote coupling 
strength Ji (J 2 ). Thick bonds have negative sign relative 
to the thin bonds, (right) Next nearest neighbor interaction 
{ 12 ) that is only allowed when the 6-fold rotation symmetry 
is broken down to 3-fold. Thick red bonds have negative sign 
relative to the thick black bonds. 


The PSG found for the vison modes can be represented 
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by these matrices. 


h = 


( ® 


T — _ 

” " 73 


,27r/6 

- i ^/2 

0 


^ 26 * 577/6 

_gi7r/6 

0 


0 0 

0 0 

g-i77/6 




V 

0 

0 

iV2 


_e-i^/e y 


0 

-1 



/ _gi77/3 

0 

0 

0\ 

1 0 

0 

0 

-1 

R3 = 

0 

1 

0 

0 

-1 

0 

0 

0 

0 

0 

gj2ir/3 

0 

\ 0 

-1 

0 

0 / 


^ 0 

0 

0 

1/ 


(D2) 


These matrices generate a group isomorphic to (C 3 x 
SL{2,’L'i)) K ( 72 ,®^ which is a finite subgroup of 0(4) 
with 144 elements. Note that this is identical to the 
group structure found in Sec. IV D. To sixth order, the 


effective Lagrangian that is invariant under this group is 
£ = 

+ a p® cos 601 — p® cos 602 — ^V^p\p 2 cos 3(0i + 02 ) 

+ 6 p® - 9piP2 + 9piP2 - p 1 + 472p®p2 cos 3(0i - 02) 

(D3) 


The fourth order Landau functional does not break 
continuous symmetries and the sixth order polynomials 
are necessary. Symmetry breaking patterns for vison con¬ 
densation are depicted in Fig. 14. When J 2 > Ji, the in¬ 
teraction prefers visons at the centers of down triangles 
to those of up triangles, and up triangles have zero weight 
in all of the soft modes. This is allowed because the hard 
core dimer constraint has been relaxed. However we can 
still study the symmetry breaking patterns. As shown 
in Fig. 14, the four distinct phases break the same sym¬ 
metries that were broken in Fig. 11 and are equivalent 
phases. 

Alternatively, we can include a further neighbor inter¬ 
action that is only allowed in the case with anisotropy. 
This is shown in Fig. 13 on the right. This interaction 
has the same symmetries as the A interactions in Eqs. (5) 
and (9). This gives the same result as above. 
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FIG. 14. (Color online) Mean-field phases diagram of the 
Landau theory [Eq. (D3)]. We have set r = —1, u = 10, and 
V = 10. The four phases are arranged to parallel Fig. 11 by 
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constraint, they tell us the symmetry breaking patterns that 
occur and is consistent with Fig. 11. 
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At the mean-field level, the dimer density deviation, d — 
d, is computed by using the expectation value, {o'fj} = 
^ap{4‘a:}{<t’i3)- Here, the mean value of the vison field, (<(»a), 
is obtained from the mean-field, Ta/f, and Eq. (17). In 
the VBSic phase, {<f)a} is zero at every other sites on the 
dual triangular lattice so that {o'fj) = 0 at every fink on 
the honeycomb lattice. Consequently, the dimer density is 
uniform on the kagome lattice. This also applies to the 
VBS 2 C phase. 





